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CFD  Validation  for  Contaminant  Transport  in 
Aircraft  Cabin  Ventilation  Flow  Fields 

INTRODUCTION 


The  risk  of  release  of  a  noxious  material  into  a  com¬ 
mercial  aircraft  cabin  during  flight  is  signiflcant  due  to 
the  worldwide  terrorist  threat  targeting  civilians.  Reac¬ 
tion  management  of  this  potential  disaster  must  be  based 
on  quantitative  analyses  of  potential  scenarios  affecting 
the  cabin  environment.  As  opposed  to  in-flight  testing, 
which  is  prohibitively  expensive  and  logistically  awkward, 
a  mathematical  modeling  with  validation  component  has 
been  proposed. 

Development  of  an  accurate  simulation  of  the  flow  held 
in  a  cabin  space  utilizing  computational  fluid  dynamics 
(CFD)  is  a  signiflcant  challenge.  First,  the  numerical 
analysis  requires  tremendous  computing  power  to  solve 
the  sets  of  equations  in  a  reasonable  time  frame.  Second, 
the  complexity  of  the  environment  must  be  handled.  Any 
obstructions  have  to  be  taken  into  consideration  since  they 
can  influence  the  flow  held  and  add  to  the  computational 
task.  For  these  reasons  initial  steps  towards  the  develop¬ 
ment  of  a  transport  aircraft  cabin  model  approach  were 
made  using  the  747  Aircraft  Environmental  Research 
Facility  (AERF,  Figure  1)  at  the  FAA  Civil  Aerospace 
Medical  Institute. 

METHODS 

The  CFD  approach  is  potentially  viable  for  predicting 
the  3-dimensional,  very  low  (2-5)  air-change  per  hour 
(ACH)  flow  flelds  associated  with  aircraft  environmental 
control  system  (ECS)  operation  in  commercial  transports. 
A  CFD  model  has  been  identifled  describing  the  unsteady, 
time-accurate,  buoyant  ventilation  flowfleld  in  an  aircraft 


Figure  1 .  Boeing  747-100  modified  for  use  as 
Aircraft  Environmental  Research  Facility  (ERF).  The 
facility  can  be  modified  to  meet  a  variety  of  research 
and  educational  requirements. 


cabin  at  cruise  conditions,  capable  of  execution  using 
an  augmented  laminar  Taylor-stabilized  flnite  turbulent 
kinetic  energy  (TKE)  model.  CFD  simulations  were  con¬ 
ducted  using  a  flnite  element  theory  implementation  of 
this  model,  via  the  University  of  Tennessee  CFD  Lab 
PICMSS  (Parallel  Interoperable  Computational  Mechan¬ 
ics  Simulation  System)  code.  Comparison  simulations 
have  used  an  available  commercial  code  {Fluent). 

The  comparative  validation  data  were  collected  utilizing 
principally  3D  sonic  anemometers  (Campbell  Scientific, 
CSAT3)  in  an  isolated  cabin  segment  of  the  747-AERF, 
which  was  devoid  of  seats.  Figure  2  provides  an  over¬ 
view  of  the  data  collection  approach.  The  selected  sonic 
anemometers  simultaneously  measured  the  3  Cartesian 
components  of  the  velocity  held  and  were  sensitive  to 
velocity  magnitude  as  low  as  0. 1  cm/sec.  The  experimental 
protocol  was  established  to  guarantee  freedom  from  drift 
over  the  lengthy  data  collection  period.  Additional  hot 
wire  data  were  taken  in  the  immediate  vicinity  of  the 
supply  jet  entry  into  the  cabin.  Flow  uni-directionality 
and  speed  in  the  immediate  vicinity  of  the  inlet  made 
this  experimental  methodology  viable. 

Data  were  collected  in  parallel  sections  that  ran  the 
length  of  the  segment  utilized.  At  each  transverse  plane 
station,  the  data  acquisition  period  was  7  minutes  with 
recording  rate  of  30  samples/second.  Within  the  3-5  min¬ 
ute  window,  the  sampled  data  were  resolved  into  mean 
and  fluctuating  components  on  1-second  intervals.  The 
mean  data  were  subsequently  averaged  over  four  30-second 
periods  to  generate  the  essential  flow  field  picture. 

RESULTS 

Figure  3  shows  the  totally  averaged  transverse  plane 
velocity  vector  held  superposed  on  an  outline  generated 
from  a  CFD  mesh.  Several  z-station  data  sets  confirm 
existence  of  the  Coanda  effect,  with  the  supply  jet  re¬ 
maining  attached  to  the  underside  of  the  luggage  bin, 
hence  clearly  penetrating  the  cabin  to  the  center  luggage 
carrier.  The  resulting  circulating  flow  field  encompasses 
the  entire  cabin  half  span,  with  a  smaller  counter-rotat¬ 
ing  flow  field  generated  in  the  upper  region  between 
luggage  carriers.  The  character  of  this  flow  pattern  does 
not  exhibit  signiflcant  z-dependence,  as  expected,  based 
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Figure  2.  Overview  of  the  movement  of  the  3D 
anemometers  utilized  for  data  collection.  Coordinate 
reference  points  in  the  segment  utilized  was  X=0  in 
the  center  of  the  cabin,  Y=0  at  the  floor  of  the  cabin, 
and  Z=0  at  the  bulkhead  that  had  been  installed  in 
the  rear  of  the  segment. 


on  the  geometric  simplicity  of  the  empty  cabin  and  its 
associated  symmetries.  The  following  comments  sum¬ 
marize  the  acquired  cabin  airflow  data: 

1.  The  velocity  held  is  deflnitely  time  varying  and 
unsteady,  with  a  periodic  unsteadiness  of  the  order 
3-4  minutes. 

2.  The  mean  velocity  extrema  within  the  Coanda  jet 
region  are  ^  1  m/s,  while  an  average  speed  in  the  cabin 
proper  is  less  than  -0.2  m/s. 

3.  Therms  data  level  holds  steady  at  "0.006-0.01  m/s, 
hence  the  average  turbulence  kinetic  energy  ^  is  - 
0.0006  m"/sh 

4.  For  average  mean  flow  magnitude  of  -0. 1 5  m/s,  the 
average  turbulence  level  is  of  the  order  -6  %. 

CFD  Predictions.  CFD  predictions  for  the  cabin  interior 
flow  fleld  were  conducted  using  2  different  code-closure 
model  systems.  The  computations  were  3  dimensional; 
however,  z-symmetry  plane  boundary  conditions  pro¬ 
duced  flow  flelds  lacking  z-dependence.  Hence,  CFD 
solutions  on  all  z-planes  are  identical,  the  transverse  plane 
discretization  of  which  is  shown  in  Figure  4  in  a  global 
view  including  magniflcations  to  illustrate  supply  and 
exhaust  ducts  meshing  details.  These  local  geometries 
were  directly  measured  for  the  CAMI  747  AERF.  The 
highly  non-uniform  transverse  mesh  contained  - 14,000 
node  points  per  plane.  The  mesh  was  identical  in  every 
z-plane. 

As  stated,  design  of  the  747  AERF  ventilation  supply 
relies  on  the  Coanda  effect  to  maintain  supply  jet  attach¬ 
ment  to  the  overhead  luggage  carrier  until  the  flow  reaches 
the  first  break  in  the  carrier  geometry.  The  experimental 
data  clearly  show  this,  and  the  supply  jet  penetrates  to 
the  center  of  the  cabin  until  encountering  the  face  of  the 
central  luggage  carrier. 

The  CFD  prediction  generated  using  the  time-accurate 
CFD  formulations  stabilized  with  a  modifledTaylor  finite 
element  approximation  (Williams  and  Baker,  1996),  as 
implemented  in  PICMSS code  (Wong  et  ah,  2002),  cap¬ 
tures  the  flow  phenomena  for  Re/L  =1000.  The  flow  fleld 
does  not  come  to  a  steady  state  after  2000s  of  simulation 
but  oscillates  about  the  essential  velocity  vector  fleld  as 
shown  in  Figure  5  a.  Figure  5  b  contains  a  magnification 
of  the  velocity  fleld  in  the  supply  duct  region  exhausting 
into  the  cabin,  and  flow  fleld  attachment  to  the  luggage 
carrier  bottom  is  essentially  steady. 

Predictions  for  comparison  were  made  using  the  Flu¬ 
ent  commercial  code  -  finite  volume  implementation  of 
a  non-time  accurate  CFD  theory.  The  Fluent  simulation 
used  the  standard  TKE  closure  model.  The  solution  at- 
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Figure  3.  The  averaged  transverse  plane  velocity  vector  field. 


Figure  4.  Cross  section  of  the  cabin  discretization  with  close  up  views  of  the  supply  and  exhaust  ducts. 
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Figure  5.  a)  PICMSS  FE  simulation  at  a  Re/L  =  1000.  b)  Velocity  vectors  in  the  area  of  the  air  inlet  supply 
region  illustrating  Coanda  effect  flow  attachment  to  the  underside  of  the  luggage  bin. 
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Figure  6.  Steady  flow  field  solution  from  the  commercially  available  code  Fluent. 
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rained  steady  state  using  the  third-order  stabilization 
scheme.  The  supply  jet  penetrates  between  the  luggage 
carriers,  Figure  6,  and  the  recirculation  region  spanned 
the  entire  cabin  region;  hence  this  flow  fleld  resembles 
the  averaged  experimental  data. 

DISCUSSION 

It  would  seem  prudent  that  protective  strategies  against 
the  threat  of  a  chemical  or  biological  weapons  attack  on  a 
commercial  transport  aircraft  should  be  developed.  Due 
to  the  complexity  of  the  cabin  environment  (airframe, 
flight  altitude,  etc.) ,  j  udgments  on  any  response  to  a  given 
situation  are  subjective  at  best.  Time  accurate  simula¬ 
tions  will  allow  a  variety  of  emergency  circumstances 
to  be  evaluated  in  a  quantitative  manner.  Insight  into 
the  best  combination  of  actions  to  take  in  face  of  a  par¬ 
ticular  challenge  can  be  built.  This  approach  also  allows 
the  response  to  any  threat  that  may  be  anticipated  to  be 
integrated  with  current  emergency  protective  equipment 
and  procedures  as  closely  as  possible.  This  would  assist 
in  making  protective  systems  aboard  aircraft  as  efficient 
as  possible  in  terms  of  both  cost  and  function. 


Accurate  simulations  are  also  applicable  to  day-to-day 
operations.  Cabin  air  quality  is  a  very  important  issue  for 
both  the  passengers  and  crewmembers.  Advanced  models 
have  the  potential  to  be  used  for  optimizing  aircraft  en¬ 
vironmental  control  system  (ECS)  parameters  for  both 
current  and  future  aircraft.  They  could  also  provide  the 
ability  to  study  the  specifics  of  a  contamination  accident 
after  it  has  occurred.  Both  quantitative  details  and  remedial 
steps  could  be  examined. 

CONCLUSIONS 

A  CFD  model  has  been  developed  to  predict  the  3- 
dimensional  commercial  aircraft  cabin  flow  field  charac¬ 
teristic  of  a  Boeing  747.  Validation  data,  as  obtained  in 
this  configuration,  confirm  the  time-accurate,  optimally 
stabilized  results  predicted  via  the  PICMSS  implementa¬ 
tion.  The  CFD  algorithm  designed  for  this  application 
produced  an  unsteady  prediction  in  agreement  with  data. 
Results  obtained  using  a  commercial  code  also  produced 
an  average  flow  field  comparable  to  the  experimental  data 
and  the  PICMSS  implementation. 
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